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We discuss the physical effects causing a modification of resonance masses, widths and even shapes 
in a dilute hadronic gas at late stages of heavy ion collisions. We quantify the conditions at which 
resonances are produced at RHIC, and found that it happens at T at 120 MeV. Although in the pp 
case the “kinematic” effects like thermal weighting of the states is sufficient, in A A we see a clear 
effect of dynamical interaction with matter, both due to a variety of s-channel resonances and due 
to t-channel scalar exchanges. The particular quantity we focus mostly on is the p meson mass, for 
which these dynamical effects lead to about -50 MeV shift, on top of about -20 MeV of a thermal 
effect: both agree well with preliminary data from STAR experiment at RHIC. We also predict a 
complete change of shape of /o(600) resonance, even by thermal effects alone. 

I. INTRODUCTION 


The first experimental data from the Relativistic 
Heavy Ion Collider (RHIC) at Brookhaven National Lab¬ 
oratory have opened a new page in studies of hot/dense 
hadronic matter. The magnitude of collective effects 
strongly suggests rapid equilibration in the system. 
Strong jet quenching, its angular dependence and the ab¬ 
sence of compensating jet all suggest the surface emission 
of jets. Why the bulk of the matter is effectively black, 
even for ~ 10 GeV hadrons/jets, remains unknown at 
this time. Another unsolved RHIC puzzle is the apparent 
large fraction of baryons at large pt . All those phenom¬ 
ena at large transverse momenta are presumably related 
to properties of extremely hot/dense hadronic matter, 
preceding formation of the Quark-Gluon Plasma. 

In this paper we would however approach the phenom¬ 
ena from the opposite side, focusing instead on the most 
dilute stages of the collision close to the freezeout. It has 
been argued over the years that under such conditions 
hadrons and especially short-lived resonances should be 
modified, with shifted mass, increased width and even 
significantly changed shape. The most dilute stage of the 
collisions, known as kinetic freezeout, is much more dilute 
than nuclear matter. In addition, most of the particles 
are Goldstone bosons, 7 r, K, 77 , which interact weakly at 
low energies. At such low density matter the proposed 
modifications are expected to be small, although quite 
observable. They should naturally be studied to set a 
benchmark, before extrapolation to denser matter close 
to the QCD phase transition, at which more species of 
hadrons are actually produced (the chemical freezeout). 

Recently the STAR experiment at RHIC reported ob¬ 
served signals of a variety of mesonic resonances such as 
p, 10 , A'*,/o(980) [1] and (f [2] as well as baryonic reso¬ 
nances (which we will not discuss in this work). A data 
sample reported corresponds to mid-central (40 to 80 per¬ 
cent of the hadronic cross section) AuAu collisions at 
y/s = 200 GeV is shown in Fig.l. (They are not yet ac¬ 
ceptance corrected, so the numbers extracted from them 
should not be directly compared to our calculations be¬ 
low.) One can see from comparison of pp and AuAu 


spectra that the shape of the p — u> pair of resonances 
in the 7 r + 7 r _ channel is quite different from its classic 
shape in e + e _ —> 7 r + 7 r _ or r decays. Unlike those lep- 
toproduction reactions, in which the masses of p and oj 
are nearly degenerate, in pp and even more so in AuAu 
collisions the p mass is shifted downward, while the to 
mass stays the same. The (preliminary) fits [1] gave 
Trip = 0.698 ± 0.013 GeV in AuAu and 0.729± 0.006 GeV 
in pp, with about unchanged width. 




FIG. 1. The invariant mass distributions in pp and 
mid-central AuAu of the 7r + 7r _ system, with a transverse mo¬ 
mentum cut 0.2 < pt < 0.9 GeV. The lines indicate contribu¬ 
tion of particular resonances according to some model we do 
not discuss in this work. 
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A. esonance mass modification: theory 

In general, resonance modification in hadronic matter 
results from mutual interaction, and is analogous to the 
well known phenomena in atomic physics known as spec¬ 
tral line shifts and broadening. 

In the first order in density of some a hadronic species 
rii , a modification of a hadron of kind j can be ex¬ 
pressed in terms of their forward scattering amplitude 
Mij(t = 0, s). It is similar to many known phenomena 
such as e.g. a modification of the photon dispersion law 
when it propagates though glass. Note that the scatter¬ 
ing amplitude is complex, and that this approach gives 
both the real and imaginary part of the dispersion law 
modification, also known as the optical potential. Note 
also that if the dependence on momenta of both parti¬ 
cles is strong, appropriate integration with the thermal 
weights should be performed. 

Unfortunately, only partial information about the scat¬ 
tering amplitude can be obtained from the experiment. 
There are two major theoretical approaches to the issue 
discussed in literature, to be called an s — channel and a 
t — channel one. 

The former approach assumes that the scattering am¬ 
plitude is dominated by s-channel resonances which are 
known to decay into the i + j channel. For most mesons 
such as 7 r, to, p , K in a gas made of pions such calcula¬ 
tion has been made [5-8] It is usually true that the sin¬ 
gle most important part of scattering amplitude is the 
lowest resonance. A general rationale for that is that 
higher resonances typically decay into states with higher 
multiplicity, with the 2-body channel in question hav¬ 
ing a small branching ratio. In the case we will discuss 
most below - the pit scattering ■ it is the axial resonance 
a-[. Generic diagrams for such processes look like that 
depicted in Fig.2(a): and in dilute matter the modifica¬ 
tion of the p is simply proportional to the appropriately 
weighted pion density (see below). 

Such an approach is of course rather limited, applica¬ 
ble only for dilute matter. Indeed, if one would like to 
follow it to second order in density with a double collision 
Fig.2(c), one should be able to extrapolate the scatter¬ 
ing data off-shell (an intermediate line with an asterisk). 
The same problem appears if one may think about loop 
corrections. Furthermore, generic 3-body collisions can¬ 
not be accessed experimentally, only their part due to 
intermediate resonances with the known 3-body decays 
(see Fig.2(d)) can. 

The predicted [5] p mass shift due to the ai in an 
equilibrated* pion gas was found to be small, not ex¬ 
ceeding -10 MeV. The sign of this effect can be explained 


‘That is, with zero pion chemical potential. The calculations 
below correspond to significant p n which of course increase 
the effect. 


as follows: states of the same quantum numbers repel 
each other if mixed, and thermally populated pn states 
have total energy mostly below the mass of ai. 



(a) (b) (c) 3 3 (d) 


FIG. 2. Forward scattering amplitude for particles of type 
1 and 2 (a) can be approximated by intermediate resonance 
3 (b). Even the double 2-body scatterings (c) include an 
off-shell particle (asterisk) which cannot be calculated with¬ 
out assumptions. An exceptional case is a single resonance 4 
which has an observable 3-body decay (d). 

Similar effects due to scattering on nucleons have been 
studied e.g. by Rapp, Wambach and collaborators, see 
[9] and references therein. There is however an impor¬ 
tant difference in sign of the effect. The major resonance 
relevant for pN scattering on nucleons, iV* (1520), is not 
only below the thermally populated states but even be¬ 
low m p + tojv: thus it pushes the p mass upward. The 
lower excitation, which is the IV*(1520)IV _1 state, is at 
the same time pushed down. Another important effect 
discussed a lot in these works is significant broadening 
of both states, leading to a complicated quasi-continuum 
spectrum in place of two well separated states. 

Let us now turn to the second t — channel approach, 
attributing the mass shifts to mutual attractive interac¬ 
tion between hadrons, such as t-channel exchange of a 
scalar isoscalar a meson, see Fig.3. We first comment on 
a significant controversy related to its association with 
the issue of depletion of the chiral condensate and chiral 
restoration phenomenon. Let us try to explain briefly 
what is its status at this point. 


j 

FIG. 3. Forward scattering amplitude for particles of type 
i and j exchanging the “radial” scalar R. 

On the one hand, a modification of the quark conden¬ 
sate can be calculated at small temperatures from general 
chiral theory [14] without a problem, the result is 

< », > (T) =< « > (0)(1 - fjl - slyj + -) (1) 

The hadronic masses appear due to chiral symmetry 
breaking, and in models like NJL it happens via the so 
called constituent quark mass which is directly propor¬ 
tional to the quark condensate. Therefore naively 
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K „ <qq>* (2) 

to p < qq > 

which is sometimes referred to in the literature as a 
“Nambu scaling”. 

Here and below in-matter quantities are marked by the 
asterisk. Brown and Rho have proposed another scaling 
( see the latest work [ 10 ] for earlier references and expla¬ 
nations) 


(3) 

m p u < qq> 

with a = 1/2 for low densities, 

The problem however is that such relations are not re¬ 
spected by low-T chiral calculations. In particular, the 
0(T 2 /f%) effect which is present in the quark conden¬ 
sate ( 1 ) as well as in of the vector and axial correlation 
functions [15], is absent in physical effects like the total 
energy of the pion gas [13] or in the mass shifts, which ap¬ 
pears only in order 0{T 4 ), see e.g. [16]. More generally, 
it is clear that a t-clrannel exchange of a sigma meson 
between a pion (with a non-derivative coupling) and any 
other hadron would violate the Goldstone theorem and 
lead to a non-zero pion mass. This contribution must 
therefore cancel out with some other diagrams, as seen 
in explicit calculations in many models. 

Predictions for next order 0(T 4 ) mass shifts are model- 
dependent and there is no theoretical consensus even 
about the sign of the contribution. The 0(T 4 ) calcu¬ 
lation performed in [16] has predicted that both p, a± 
masses decrease slightly, in equilibrium T = 150 MeV 
was found to be dm p 4 = —6AMeV. 

In contrast to that, Pisarski [17] has built a set of 
simple effective models incorporating 7 r, p, a± interactions 
(which all of course respect chiral symmetry) and con¬ 
cluded that no prediction can be really made. In one of 
them, enforcing VDM at any T, the result is that the 
p mass should grow with T significantly, joining ai at 
at chiral restoration point at the following intermediate 
mass 

m 2 p (T c ) = m 2 ai m) = ^m 2 (0) + ^m 2 ai (0) = (962 MeV) 2 (4) 

In contrast to scalar exchanges with pions, such ex¬ 
changes between other hadrons definitely exist. (For nu¬ 
cleons those are in fact responsible for a large portion of 
nuclear forces and binding.) The effect of nucleon den¬ 
sity on quark condensate related with the nucleon sigma 
term [18] is large even at nuclear matter density 


< qq >* _ 1 _ S jjv n 
<qq> f 2 ml 


1 - .36- 


n 


n 0 


(5) 


where n o is the nuclear matter density. The resulting BR 
scaling prediction 

m* = mp{ 1 — 0.36— + ...) 1 / 2 ( 6 ) 

” 71 r. 


leads to good agreement with experimental information 
on scattering on nuclei. 

This apparent contradiction is discussed in multiple 
papers. We would follow [11] where it has been especially 
clearly emphasized that one may think of two different 
scalar fields which can be exchanged. One is the old a of 
the linear sigma-model, a chiral partner of the pion. The 
other we will call R t 

R 2 = a 2 + 7T 2 (7) 

The main difference between VEVs < a > and < R > 
is that the former gets the pion loop 0(T 2 ) correction 
while the latter does not since in the excitations along 
the chiral circle the radius does not change. In a different 
language, the interaction Lagrangian describing how R 
interacts with pions is different from a sigma model one, 
it has a derivative coupling like for any other particle. 
Therefore R,-exclianges with pions does not give them a 
mass and does not violate the Goldstone theorem; there 
is no need for cancellations. The mass shifts for N and 
p can then in the first order be written as R-exchange 
diagram or [ 11 ] 


SM n = 'Um p = - 9 \ ( 8 ) 

2 m R 

In this simple picture go = M n /f n ss 10 and R can be 
viewed as the Walecka scalar meson, see also [12] where 
R is denoted as %• 


B. Resonance mass modification: experiments 

We will not go into long history of detailed p-meson 
shape, but simply state that the difference between its 
appearance in elementary reactions (e + e~ —> 7 r + 7 r _ or 
r —y v T ir + ir~ decays) and hadroproduction reactions is 
by now well known. The latest Review of Particle Prop¬ 
erties (RPP) [19] averages the rho mass for these two 
classes of experiments separately, with a clear systematic 
difference of the order of 10 MeV between them: 

m leptoproduced = 775 g ± 5 Me y^ ( 9 ) 

m hadroproduced = jqq 5 ± Me y ( 10 ) 

. We return to discussion of this effect in section IIB 
below. 

A p mass shift Sm p = 160 ± 35 MeV observed in nu¬ 
clear experiments has been reported by Lolos et al [3] in 


^From the radius - in [11] it was called 0 which may be 
mistakenly taken as an angle or may be confused with the 
widely used 0 field, the U(l) pseudoscalar. R is the x of 
Brown and Rho [12]. 
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the 3 He( 7 , p°)ppn reaction. An even larger shift has been 
reported by Huber et al [4] when even lower 7 energies 
were used. Although 3 He is not a dense nucleus, with 
average density of about half of nuclear matter, the kine¬ 
matic of this experiment is sub-threshold for one nucleon. 
It was therefore argued that the emerged p is actually in 
the vicinity of at least two nucleons, so that the local 
density is large. 

Heavy ion dilepton experiments at SPS, such as HE- 
LIOS3,NA38 and especially CERES have indeed found 
that the vector spectral density of hot/dense matter is 
very different from that observed in elementary pp col¬ 
lisions. While in the latter case there exist a clear gap 
between p, w mesons and low-mass pairs originated from 
Dalitz decays of 77 and other light mesons, no such gap is 
present in heavy ion collisions and a continuum of excited 
states is seen down to the invariant mass of ~ 400 MeV. 
Although these results agree with various calculations 
such as the Rapp-Wambach ones mentioned above, these 
dilepton data are still very limited statistically and have 
very crude mass resolution which does not make it pos¬ 
sible to see the p modification or even separate it from 
the u! or even the N*(1520)N~ 1 . 

The situation is very different for STAR experiment 
we mentioned in the Introduction, in which the shifts are 
smaller and the width apparently is not growing, while 
the invariant mass resolution is excellent and all reso¬ 
nances are seen as clearly resolved separate peaks. 


II. MODIFICATION OF p MESON AT RHIC 
A. Kinetic freezeout at RHIC 

Early papers on statistical and hydrodynamical models 
of multiparticle production had assumed complete equi¬ 
librium matter untill the final freezeout. Clear separa¬ 
tion between two different stages, the chemical (at T c h) 
and kinetic (at T k ) freezeouts has been made in the mid- 
1990’s. Basically it reflects very different rates of elastic 
and inelastic reactions at low energies dominating rescat¬ 
tering of secondaries. 

In particular, in ref. [21] the following argument has 
been put forward: at AGS/SPS energies most of the ob¬ 
served velocity of collective radial flow can only develop 
in between the two, at T ck < T <T k . 

In that paper also the following natural definition of 
the kinetic freezeout condition has been given: a sec¬ 
ondary (say a pion) emitted at this time has equal 
chances to be either rescattered or escape: 

Pescape — Prescattering — ^ (H) 

which we will be using below. 

At the particular conditions we are interested in, at 
RHIC t/s = 200 GeV, the numerical values of the pa¬ 
rameters are as follows. The rapidity density dN/dy at 


mid-rapidity of central AuAu collisions for relevant par¬ 
ticles are about 300 for each species of pions, 29 protons 
and 22 antiprotons. These numbers have large feed-down 
from resonance decays, which should be accounted for 
when their ratios are used to derive the chemical freeze¬ 
out conditions. From [1] it was determined that those 
are about 


T ch = 170 MeV p b = 28 MeV (12) 

By definition of chemical equilibrium, all other chemi¬ 
cal potentials for non-conserved charges are zero at this 
stage. 

After inelastic reactions are frozen and only elastic re¬ 
scatterings take place, the number of particles of each 
species are separately conserved. Therefore chemical 
potentials for all species /j, become non — zero at the 
hadronic stage, see details in refs [22,21,23,24]. The only 
important numbers for our discussion below are the value 
of the kinetic freezeout temperature and the values of 
chemical potentials for pions and nucleons. Assuming 
entropy conservation (adiabatic expansion) and natural 
conservation of baryons one can determine all parame¬ 
ters along a particular cooling path, provided one point 
on it is known. Using (12) one gets it, and using either 
kinetics, or hydro-based fit to data, one finds that the 
kinetic freezeout at RHIC happens at 


T k « 100 MeV, p n « 81 MeV 
p N « 380 MeV p K « 167 MeV (13) 


This translates into the following densities of the pions 
(all 3 of them) 

n w «0.06/m -3 (14) 


With p/p « 0.75 we get the following density of all nu¬ 
cleons plus anti-nucleons together 

Ujv+jv ~ 0.0075 /to -3 (15) 


or only 1/20 of the nuclear matter density. Accounting 
for all baryonic resonances will increase them by nearly 
a factor 3 sand bring them in agreement with measured 
final yields mentioned above. 

The next issue we discuss is how the freezeout condi¬ 
tions for most particles (pions) relate to those for reso¬ 
nances. When those are actually observed as peaks in 
the invariant mass distribution of several secondaries, it 
is only possible when none of the secondaries have been 
rescattered. Suppose for simplicity the resonance is pro¬ 
duced and decays into the same mode with N identical 
secondaries of type i (N=2 for p —> 2n). Its observ¬ 
able production at time t is proportional to the following 
combination of the density needed for production and the 
probability that all secondaries do escape, namely 


P{t) 


n 2 


1 N 

pOO 

( t)exp (— / dt' E nj(t') < aij V ij » 

Jt + T A 


(16) 
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where n,(i) is time-dependent density of relevant sec¬ 
ondary in matter, and < (JijVij > are thermally averaged 
cross section and velocity for scattering on type-j particle 
and t is the resonance lifetime. 

This expression can be simplified significantly, pro¬ 
vided the lifetime is ignored, and the matter expansion 
is parametrized by a single power behavior rii(t) ~ l/f a . 
The cross section is approximately time-independent and 
can be eliminated in favor of the freezeout time, using its 
definition (11), leading to 


P, 


resonance 


(t) 


t a exp(—ln(2){—) 


tk , a _i 


N 


(17) 


which for any N has the maximum at 

Ues = tk [ln(2) * (1 - l/a)]“ (18) 

As hydro simulations show^ the value of the index is actu¬ 
ally not constant over time; it changes from the Bjorken 
a=l at early time of Id expansion to a value a ss 2.5 at 
relevant times t ~ 10 fm/c , not quite reaching ultimate 
3 for 3d expansion. For this index, t res « 0.61*,. The dis¬ 
tribution (17) for N=2 is plotted in Fig,4, and one can 
see that the peak is rather sharp. 



t 

FIG. 4. Probability distribution over the 2-pion resonance 
production time t, in units of pion freezeout time tk- 

Finally, one can deduce from this time the relevant 
temperature using a power-fit for density at the hadronic 
stage n ~ T 5 which gives 


T r , 


T k 


tk 


—a /5 


1.2 


(19) 


In effect, we found that the observed resonances are pro¬ 
duced somewhat earlier and a bit higher temperature 
T res compared to the kinetic production time, for all N. 
The conditions on the same adiabatic path at such tem¬ 
perature are: 


T res « 120 MeV, p n « 62 MeV 

p N « 270 MeV p K w 115 MeV (20) 

with densities about 1.4 higher than at the kinetic freeze¬ 
out. 


B. Elementary effects due to the heat bath 

Processes like 7T7 t scattering in an equilibrium heat 
bath can be represented as follows 

dW = [f(pi)f(pi)\ \M\ 2 dT [(1 ± /(p 3 ))( 1 ± /(p 4 ))] (21) 

where (i) the first two factors are thermal occupation 
factors for the incoming particles; (ii) M is the matrix 
element of the re-scattering, (iii) is the phase space and 
finally (iv) represent final state corrections for induced 
radiation, ± for final bosons/fermions; (iii) In dilute mat¬ 
ter one can presumably use unmodified matrix element 
and the phase space. 

The socalled “elementary” hadronic reactions like pp 
collisions at high energies are in fact very complex phe¬ 
nomena, which are far from being understood. Unlike 
heavy ion collisions, those do not show collective phe¬ 
nomena like flow and one may think they certainly do 
not contain equilibrated hadronic matter. And yet, as 
noticed back in 1960’s by Hagedorn and others, the par¬ 
ticle composition in pp and even e + e— collisions can 
be well reproduced by statistical models, see e.g. recent 
comparison [20]. The fitted temperature is about the 
same as the chemical freezeout temperature T c h men¬ 
tioned above. Two- and even many-body distributions 
follow thermal distributions very well, as data showed al¬ 
ready in 1970’s. So it is natural^ to apply it to 2 pions 
which are forming the p meson in question as well. 

Now, if two initial pions come from the heat bath even 
in pp, the initial thermal factors should be present. At 
the time the resonance decays, on the other hand, the 
secondaries making that heat bath are already sufficiently 
far away. From this it follows that (i) the resonance decay 
happens in vacuum, its matrix element is unmodified, 
and (ii) the two final state induced radiation factors are 
absent. 


*We thank P.Kolb who provided plots from which the index 
value was estimated. 


§The authors thank P.Braun-Munzinger who suggested this 
idea to them. 
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The Boltzmann factor is of course a function of the en¬ 
ergy in the matter rest frame, not in the resonance rest 
frame. However as a simple first approximation one may 
neglect the thermal motion of rather heavy resonances 
and assume them to be at rest in the heat bath. This 
simply combines the Boltzmann factor exp(—M/T) con¬ 
taining the invariant mass M 2 = (p_|_+p_) 2 with the (ap¬ 
propriately modified for p-waves if needed) Breit-Wigner 
matrix element. 

The position of the maximum can then be obtained in 
a simple analytic way, provided the resonance width is 
large enough. Approximating the resonance by 

exp(—M/T) _ 

1 + (4/T 2 )(M — TTlres) 2 ~ 1 ’ 

a o MV 2 

expH^XM 2 - 2 Mm res - M 2 + 

one gets the mass shift due to (sufficiently large) temper¬ 
ature T to be approximately 


a 


31 


° O Jt 

° o 

<b 
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6m 


max 


r 2 

8 T 


(23) 


16l 


For p emitted at T = 120 MeV this expression gives 
the shift - 23 MeV. This is to be compared to -10 MeV in 
average low energy hadroproduction reactions (10) and 
to a shift of about -30 MeV in preliminary STAR pp data 
[1] mention in the Introduction. 

The important point here is that this effect is quadratic 
in width. For example, the 3 times more narrow K* is 
predicted to be shifted by this effect an order of mag¬ 
nitude less. In contrast to that, much wider resonances 
- e.g. the famous cr-meson, now listed as /o(600) [19] 
with a width of about 300 MeV - are changed beyond 
recognition. Applying the same expression we see that 
in the pp case a Boltzmann factor at T ~ 160 MeV trans¬ 
form it into a very wide structure (which is indeed seen 
in hadroproduction), while at kinetic freezeout of heavy 
ion conditions T ~ 100 MeV it peaks at very low masses 
instead, see Fig,5(b). (Since it is an s-wave, not p-wave 
resonance as p, its threshold suppression is weaker.) 

We argued above that in the pp case the 2 pions form¬ 
ing the p come from a heat bath, while the outgoing 
pions propagate essentially in free space. In contrast to 
that, in heavy ion collisions the local expansion rates are 
slow compared to the p lifetime, and therefore it decays in 
matter. So, in principle one should include the final state 
factor (1 ± /) 2 . Direct experimental access to the magni¬ 
tude of thermal occupation rate at freezeout is provided 
by the integrated HBT correlator. For p it is however a 
small effect modifying the width by about 10 percent and 
not affecting the mass. It is larger for u>: see discussion 
of it in [5]. 
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FIG. 5. The shape of the p (a) and /o(600) (b) as they ap¬ 
pear in elementary reactions (squares), modified by a Boltz¬ 
mann factor at chemical freezeout T=165 MeV (diamonds) 
and at chemical freezeout (circles) 


C. Modification due to rescattering on pions, kaons 
and nucleons 


The formulae are all well known [5], for the record we 
use the resonance 2-pion formfactors in the vector and 
scalar channels, which include the phase space and p 2 in 
the p-wave matrix element in the width 


FF sq 


m^ 

(m 2 — E 2 + f) 2 + m 2 r v 2 


(24) 
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FF„, 


m* 


where 


(m 2 - E 2 + f) 2 + m 2 T a 2 
r mk 3 


r„ := 


r„ := 


Ekv 3 

r mk 
E kv 


k := — \/E 2 — 4 m n 


(25) 

(26) 

(27) 

(28) 


kv := \Jm 2 — 4 m K 2 (29) 

We found that the resonance mass shifts due to s- 
channel resonances depend rather weakly on the emis¬ 
sion time, making predictions rather stable. In particu¬ 
lar, the difference between (20) and (13) is at the level 
of 10-15 percent, much less than density difference which 
is 40 percent. This is explained by a specific shape of 
the real part, with negative and positive contributions 
weighted with a thermal weight. At lower temperature 
higher states contribute less and their cancellation is less 
pronounced. 

As seen from the Table, the interaction of p with the 
pions lead to negative shift of the order of -30 MeV, while 
with kaons it is basically zero. 

Due to IV* (1520) the interaction with nucleons leads 
to positive shift. There are however resonances above 
the threshold which are not included; those will some¬ 
what reduced our estimate. On the other hand, due to 
S=I=3/2 of the A state, it has rather large statistical 
factor and significant population at the relevant stage. 
There are sub-threshold p — A resonances in the 2 GeV 
mass region and thus one can expect an effect similar 
in sign and magnitude to those from A* (1520), further 
increasing the upward shift. 

An opportunity to test the role of Np interacton is 
provided by heavu ion collisions at lower collision ener¬ 
gies, of the order \Js few GeV/N (like expected in new 
GSI proposed heavy ion project). Since tha increases the 
baryon fraction by an order of magnitude, the combined 
effect of p — N 1 p — A resonances can compensate or even 
overcome the negative shifts discussed above. The mea¬ 
surement of the exact mass and shift of p in rnr mode at 
such conditions is therefore of significant interest. 


D. The effect of mutual attraction 

As explained above, the majority of the particles in 
the matter are Goldstone bosons 7r, K 1 p which do not 
produce attraction due to R-exchanges. There should 
however be an effect due to other particles, such as vec¬ 
tor mesons and baryons. Unlike the s-channel resonances 
discussed above, the t-channel attraction scales as den¬ 
sity, and thus is 40 percent stronger at (20) as compared 
to (13). Using a simple expression for the mass shift (8) 
and the R.-scalar mass Mf j ss 800 MeV one gets 

Snip « -28 MeV (30) 

(in agreement with that from the BR scaling) due to all 
B + B. An additional shift of the magnitude of 

5m 1 ’ « —10 MeV (31) 

comes from scalar R-exchanges between p and vector 
mesons p,u,K*. 


E. Modification of the rho width 

Finally let us briefly address the issue of the width 
modification. The main difference between the two mech¬ 
anisms of the mass shift discussed above is that the t- 
channel attraction has no imaginary part and is not as¬ 
sociated with the broadening, while the s-channel reso¬ 
nances lead to well defined broadening, of the order of 
the sum of the absolute values of the mass shifts which 
is about 50 MeV. 

On the other hand, there is a “kinematic” effect work¬ 
ing to the opposite direction. The negative mass shift 
discussed above automatically reduces the width, both 
because of the reduced phase space and also due to the 
power of p in the P-wave matrix element. The magnitude 
of this effect for the predicted mass shift is 

&7TI 

<5T P = 3—^T p « -50 MeV (32) 

m p 

So, inside the accuracy one may claim at the moment, 
these two effects seem to cancel each other. 


III. DISCUSSION AND OUTLOOK 


TABLE I. A set of resonances considered 


Name/Mass 

Width 

Branching 

Mass Shift 

oi(1260) 

400 

0.6 

-19 

02(1320) 

104 

0.7 

-15 

A'i(1270) 

90 

0.4 

+1.6 

A' 2 (1430) 

100 

0.087 

-0.4 

N* (1520) 

370 

0.2 
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Summarizing the paper, we have found that: (i) a sim¬ 
ple initial state occupation effect is sufficient to explain 
the p mass shift between pp and leptoproduction; (ii) 
the effect is predicted to be quadratic in width, so it 
should be an order of magnitude smaller for K*. At the 
same time it leads to very significant shape deformation 
for the sigma /o(600) pushing it downward toward the 
threshold; (iii) we found that all resonances are emitted 
at some time between chemical and kinetic freezeout, at 
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T w 120 MeV; (iv) dynamical effects due to interaction 
of the p with surrounding matter lead to about -50 MeV 
additional p mass shift, relative to a shift in pp mentioned 
above. 

All these conclusions about the mass shift of the p 
agree with the first STAR data shown at Quark Matter 
02 [1], Hopefully quantitative results, with fitted values 
of mass and width versus centrality, will be soon available 
for more quantitative comparison. 

The main uncertainties of our calculations of s-channel 
resonances come from poorly known branching ratios for 
modes containing the p mesons we focus on. It is clear 
that smaller contributions of many more resonances are 
missing, which may or may not average out to zero. We 
think the accuracy of our predictions are therefore at the 
level of 20-30 percent. Observation of similar phenom¬ 
ena with several resonances and at variable conditions 
(such as N/n ratio) would help to reduce this uncer¬ 
tainty. Very interesting would be precise measurements 
of the line shape of two more vector mesons, K* and <j>, 
and possibly a shape of lo in 3 pion channel. 

In this paper we have not addressed in detail the issue 
of the p width and its broadening: however as data would 
become more accurate, it would provide an opportunity 
to better test the relative role of s-channel and t-channel 
conntributions to the mass shift. 

STAR data also display a peak in 7T7 t channel below 
p , at invariant mass M ~ 600 MeV. It can be partly 
a = /o(600), which is however predicted in this work 
to be deformed and move downward. It is likely to be 
2/3 of the u> with the 7r° missing, or a trace of baryonic 
resonance like AT*(1520). 
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